1 


Stirling Analysis Comparison of 
Commercial VS. High-Order 

Methods 

Rodger W. Dyson, Scott D. Wilson 
Roy C. Tew, and Rikako Demko 
Thermal Energy Conversion Branch 

August 15, 2005 






Stirling Simulation Numerical Error 
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•Remeshing interpolation 
•Layering interpolation 
•Diffusive time advance 
•Low grid quality/skewness 
•Sliding interface interpolation 
•Artificial entropy 
•Turbulence transition 
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Schematic with Springs and Dampers 


Effect of Heating & Cooling / Artificial Transient Heat Loss 
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Dual Opposed Convertors 
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• High Efficiency - Low Mass Space Power 



• Free Piston Geometry is Essentially Smooth 








Whole Engine Simulation 
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•Bounded by Walls - No need for 
nonreflecting B.C. 

•Kolmogorov scales fairly large 

•Steep thermal gradients 

•No shocks/ subsonic/transitioning 

•High-order friendly 
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Curvilinear Features 
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Stability Analysis 
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• Courant (CFL) number, 
r = c At / Ax 

• Von Neumann number, 
v = |a At / Ax 2 

• Linear Viscous Burger’s 
Equation 




Compact Scheme 6 th Order in Space 



a=2/l \a=\21\ 16=3/1 1 
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Runge-Kutta 4 th Order 

R(u) = -cu x +juu xx 

u m =u"+ — R ” 
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W (3) = u n + A£/? 2 

w ( " +1) = u n + — (iT + 2i? (1) + 2i? (2) + R (3) ) 
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i? (1) = i?(w (1) ),i? (2) =i?(w (2) ),i? (3) = i?(w (3) ) 
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Current Practice 
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• Implicit 1 st or 2 nd order in time commercially 

• 1 st or 2 nd order in space implicit 

• Explicit/implicit 4 th order in time academically 

• Implicit 6 th order compact scheme in space 
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CFL 
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Compact Scheme Stability Range 

• Domain size affects stability since implicit 






Basic UHF Technique 
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Derivatives of Cell Averages 
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Grid Singularity Resolution 
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C4o0 Linear Viscous Burger’s Equation 
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Efficiency Improves with Accuracy 
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Computational Domain Schematic 
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Amplification Factor Comparison 
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Overview of Segregated Solution 
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Heat Transfer Test 



Reduces to: 
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Temperature 


Comparison of Commercial & Advanced 
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Exact Solution 
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Commercial Comparison 
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Turbulence Transition Efficiency 
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Points per Kolmogorov Wavelength 
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Wave Equation Amplification 
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Aliased Frequency Amplification 
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Conclusions 
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• Low Reynold’s number, wall bounded flow 
allows economical use of large eddy 
simulation for turbulent transition modeling 

• UHF and Compact comparable at conjugate 
heat transfer 

• UHF much better for turbulence modeling 

• Modern methods much more efficient than 
those currently available commercially 
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